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Abstract 

We investigate CP-violation effect in the long-baseline neutrino oscillation in 
the four-neutrino model with mass scheme of the two nearly degenerate pairs sep- 
arated with the order of 1 eV, by using the data from the solar neutrino deficit, the 
atmospheric neutrino anomaly and the LSND experiments along with the other 
accelerator and reactor experiments. By use of the most general parametrization 
of the mixing matrix with six angles and six phases, we show that the genuine 
CP-violation effect could attain as large as 0.3 for AP(z/^ 1/^) = Pii^^ —>■ 
Ur) — P{i7fi i7r) and that the matter effect is negligibly small such as at most 
0.01 for AP(iy^ i/^) for Am^ = (1 — 5) x 10^^ eV^, which is the mass-squared 
difference relevant to the long-baseline oscillation. 
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I Introduction 

It has long been assumed that neutrinos are massless. However, since the atmo- 
spheric neutrino anomaly was discovered by several experimental Collaborations IH] 
and was affirmatively confirmed by the Super-Kamiokande[0, people have come 
to think through the neutrino oscillation interpretation for the anomaly that neu- 
trinos seem to have a certain amount of mass. Together with the solar neutrino 
deficit 0, the anomaly has been analyzed in the three- neutrino model |^ and two 
typical mass scales have been derived for the neutrino mass-squared difference; 
^^atm — (0-5 — 6) X 10"'^ eV^ with a large mixing angle of sin^2^^atm > 0.82 as 
the z/^ Uj. oscillation from the atmospheric neutrino anomaly and Aml^i^^ = 
(10^^^ — 10"^) eV^, a large range depending on the three solutions of the vac- 
uum oscillation and the MSW solutions in the matter with small- and large-angle 
mixings from the solar neutrino deficit . 

As in the quark sector, CP violation would be a characteristic feature in the 
three-neutrino model. It has been shownQ[^ @@ by using the constraints on 
the mixing matrix elements obtained from the analyses of these anomalies along 
with the results from the other accelerator and reactor experiments that the CP 
violation effect, defined as a difference of the oscillation probabilities between 
the neutrino and the antineutrino, is typically 1 — 3% even in the long-baseline 
neutrino oscillations, depending on the assumed mass hierarchies. 

On the other hand, sterile neutrinos were considered in the context of neutrino 
oscillations[0][jlT|. After that, a four-neutrino model of the ordinary three active 



neutrinos and one sterile neutrino was introduced with a mass pattern of two 
nearly degenerate pairs separated with a mass gap of the order of leV motivated 
from the hot dark matter |12[ and then, by using the only one possible positive 
evidence from the terrestrial LSND experiments on the oscillations z/^ Ve and 
~^ [01' which suggest the mass scale of AmLSNo = (0-3 — 2.2) eV^, the 
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four- neutrino model with the same mass pattern as the above is studied [14-20]. 
In this model, a sizable CP violation effect is shown to be possible in the long- 



baseline experiment s ||15|| |2l|l , and different magnitudes of the probability difference 
between the CP-conjugate channels are expected in between the three- neutrino 
model and the four-neutrino model |^ by using the most general parametrization 
of the mixing matrix[^]. And, some features of CP asymmetry defined as the 
normalized probability difference are discussed in the long-baseline experiments 
at a neutrino factory |23| . 

We will investigate the CP violation effect in the long-baseline neutrino oscil- 
lations numerically in more detail in the four-neutrino model with mass scheme 
of the two nearly degenerate pairs separated with the order of 1 eV by using the 
most general parametrization of the mixing matrix, and in addition we will study 
the matter effect in the four-neutrino model. 

The paper is organized as follows. In Sec. II the four-neutrino model we use 
here is presented and the expressions of the difference of oscillation probabilities 
between the CP-conjugate channels are given both in the exact form and in the 
approximate forms relevant to the short-baseline and the long-baseline neutrino 
oscillations for the neutrino mass scheme mentioned above. In Sect. Ill con- 
straints on the neutrino mixing matrix are derived by using the solar neutrino 
deficit, atmospheric neutrino anomaly, Bugey reactor experiment, CHOOZ exper- 
iment, LSND experiments, CHORUS and NOMAD experiments and the other 
accelerator and reactor experiments. In Sect. IV the most general parametriza- 
tion of the mixing matrix is adopted to obtain the constraints on the mixing angles 
and phases from the ones on the mixing matrix derived in Sect. III. And then, 
CP-violation in the long-baseline neutrino oscillation is investigated on the basis 
of these constraints. The behavior of the oscillation probability differences is an- 
alyzed in detail with respect to the two relevant phases of the mixing matrix and 
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/\m?/E. The matter effect is shown to be neghgibly small in the four-neutrino 
model with the mass scheme adopted here. Finally, Sect. V is devoted to the 
conclusion. 



II The four-neutrino model 

In order to consider the solar neutrino deficit, the atmospheric neutrino anomaly 
and the LSND experiment, we will take the four-neutrino model with the three 
ordinary active neutrinos and one sterile neutrino with three different scales of 
the neutrino mass-squared difference, Amg^jg^j, = (10^^^ — 10^^) eV^,Am^^jjj = 
(10-3 - 10-2) eV^ and Am^gj^o = (0.3 - 10) eVl 

Under the neutrino oscillation hypothesis [p4| p5| , the flavor eigenstates of neu- 



trinos are the mixtures of mass eigenstates with masses mj(i = 1, 2, 3, 4) as follows, 

4 

T^a = ^UaiVi, a = e,fi,T,S (1) 



1=1 



where z/g, t'^ and Ur are the ordinary neutrinos and Ug is the sterile neutrino, and 
U is the unitary mixing matrix. The neutrino oscillation probability of z/q, 
in vacuum is given in the usual manner in the four-neutrino model by 

Pii^a ^ i^p) = 5a(3-4: J2 MKkUajU^jU^k) sin^ Afc,+2 J2 MKkUajU;^Upk) sin2A 

k>j k>j 

(2) 

where A^j = Am\jL/{AE), L being the distance from the neutrino source and E 
the energy of neutrino. The oscillation probability for the antineutrinos is given 
by the exchange oi U ^ U* in Eq.(^. And, the probability difference between 
CP-conjugate channels given by 

= AY^lm{U:,U^,U;pp,)sm2A,, (3) 

k>j 

is a direct measure of the genuine CP-violation effect in the neutrino oscillation 
in vacuum fEBtl. 



The four neutrino masses should be devided into two pairs of close masses 
separated by a gap of about leV in order to accomodate with the solar and atmo- 
spheric neutrino deficits and the LSND experiments along with the other results 
from the accelerator and reactor experiments on the neutrino oscillation. There 
are the following two schemes for that mass pattern; (i) Am'^^Ya.r = ^"^21 ^ 
^"^atm = ^^Is < ^"^LSND = ^"^32. and (11) Am^^i^, = Am^g < Aml,^ = 
Am2i ^ Aml^^Y) = Am|2, where Am\j = m\ — m^. We will adopt the first 
scheme in the following analyses, and the second scheme can be attained only 
through the exchange of indices (1,2) (3,4) in the following various expres- 
sions such as the oscillation probabilities. In the first scheme, the measure of 
CP violation in the neutrino oscillation in vacuum is given for the short-baseline 
experiment ( LjE ^ 1 [km/GeV] ) as follows, 

AP^p ~ 4[Im(f/:3[/«2t/|2[/^3)+Im(f/:3f/«i[/^i[/^3) 

+ \^{lJl^U^2Ul2Ufi^) + \^{lJl^U.xUl^U(i^)\ sin2A32, (4) 

since A21 and A43 <^ 1, and A41, A42, A31, A32 — 1. AP^^ in Eq.(|D is zero due 
to the unitarity of the mixing matrix U . So, CP violation is negligibly small in 
the short-baseline oscillation experiments in the four-neutrino model. 

On the other hand, for the long-baseline experiment { L/E = 100 — 1000 
[km/GeV]) the probability difference in vacuum is given as follows, 

AP^p ~ 4Im([/*4[/„3t/^3^/34) sin 2 A43, (5) 

since A21 ^ 1, A41, A42, A31, A32 ^ 1, and A43 ~ 1. There are six APa/3's, 
that is, AP^e, APer) AP^^, APes; AP^s, and AP^g. Three of these six AP^/j's are 
independent due to the unitarity of U for the approximate expression of APap in 
Eq.(l^) as well as for the exact expression in Eq.(|D. 
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Ill Constraints on the mixing matrix U 

In order to numerically calculate the oscillation probability differences AP^^, 
we will derive the constraints on the mixing matrix U from the solar neutrino 
deficit, atmospheric neutrino anomaly, LSND experiments and the other terrestrial 
oscillation experiments using the accelerators and reactors. 

(i) Solar neutrino deficit 

Since A21 ~ 1 and all the other five Ak/s are enormously larger than 1, the 
survival probability of Uf. is given from Eq.(|D by 

PsoUl'e^l^e) ^ l-4|f/einf/e2|'sin2A2l-2|t/e3p(l-|t/e3|'-|t/e4r) 

-2|f/e4p(l-|f/e4|'), (6) 

where the unitarity of U is used. For the solar neutrino deficit, there are three 
different kinds of solutions, that is, the vacuum solution and the MSW solutions 
with small and large angle mixings, and a unique solution is not yet found, so 
that we will not use this deficit in order to obtain the constraints. 

(ii) Atmospheric neutrino anomaly 

Since A21 ^ 1, A43 ~ 1 and A41, A42, A31, A32 ^ 1, the survival probability of 
is given by 

(7) 

By using the data from the Super-Kamiokande experiments, that is, sin^ 25atm > 
0.82 for 5 X 10""^ < Aml^^^ < 6 x 10"^ eV^, and expecting from this data that 
\U^i\'^ + |^|j2p ^ I5 the following constraint is obtained, 

\U^3\^\U^4\' > 0.205. (8) 

(iii) The Bugey experiment^ (including Krasnoyarskg, CDHS|2| and CCFR[g 
experiments) 
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By being typically represented by the Bugey reactor experiment with L/ E = 3—20 
[m/MeV or km/GeV], since A21 < 1,A43 < 1 and A41, A42, A31, A32 ~ 1, the 
survival probability of z7g is given by 

PBu^ey{^e ^e) ^ 1 - 4(1^/^3^ + \U,^\''){1 - |t/e3|' - \U,^\'') sin^ A32. (9) 

If we use the data from the Bugey experiment conservatively, that is, sin^ 2^Bugey < 
0.1 for 0.1 < Am^ < 1 eV^, the following constraint is obtained, 

If/esT + \Uei? < 0.025. (10) 
(iv) The CHOOZ experiment |]3T|] 

This experiment is the first long-baseline reactor experiment, since L ~ 1 km and 
E ~ 3 MeV so that L/E ~ 300 [km/GeV]. Therefore, A21 < 1,A43 ~ 1 and 
A41, A42, A31, A32 ^ 1, and the survival probability of 14 is given by 

PcnOOz{l^e ^ l^e) ^ 1 -4|f/e3n?7e4r siu^ A43 -2(|f/e3|' + |f/e4r) (1 - |?7e3|' - |f/e4|'). 

(11) 

By using the data from the CHOOZ experiment, that is, sin^ 2^chooz < 0.12 for 
3 X 10"^ < Am^ < 1.0 X 10~^ eV^ and adopting Eg . (|T0|) , the following constraint 
is obtained, 

^UMUeA.? < 0.12. (12) 
If we use, however, the constraint of Eg . ([T0|) and the uneguality of 2 1 f/e3 1 1 t/e4 1 < 



|f^e3p + |f^e4p, & Constraint 4|?7e3p|t/e4p < 6.3 x 10""^ is obtained so that Eg.(|T2]) 
is included in the constraint from the Bugey experiment, 
(v) The LSND experiments [p!^ 

This experiment is of the short baseline, L/E = 0.5 — 1 [m/MeV]. Since A21 ^ 
1, A43 <^ 1 and A41, A42, A31, A32 ~ 1, the oscillation probability of z/g is 

expressed as follows, 

Plsnd(/^m ^e) ^ -4Re [{U;.,Uez + U;,Ue4){U^iU:, + U^2U:2)] sin' A32 

= 4|f/;3f/e3 + f/;4f/e4rsin2A32, (13) 
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where the unitarity of U is used. By using the data from the LSND experiments, 
that is, sin2 2^LSND = 1-5 x lO^^ - 1.0 x lO^^ for 0.3 < Am^gNo < 2.2 eV^ the 
following constraint is obtained, 

|f^;3^e3 + U;^U,^\ = 0.02 - 0.16. (14) 



vi) The CHORUS m and NOMAD |33| experiments 



These experiments are also the short baseline ones searching for the z/^ oscil- 
lation, L/E = 0.02-0.03 [km/GeV]. Since A21 < 1, A43 < 1 and A41, A42, A31, A32 
10^^ — 10^^, the oscillation probability is given by 

^CHORUS/NOMAD (Z^M ^ l^r) ^ 4:\U*.;^Ur3 + U*^Ur4\'^ slu^ A32. (15) 

By using the data from the latest NOMAD experiment, sin^2^NOMAD < 0.3 for 
Am^ < 2.2 eV^, the following constraint is obtained, 

|f^;3f^r3 + f/;4f/r4| <0.28. (16) 

Among the above-mentioned six typical phenomena and experiments, the use- 
ful constraints are of Eqs. (D, (|lOD, O and ([T|). 



IV. CP violation in the neutrino oscillations 

In this section, by using the constraints obtained in the previous section, we 
will numerically investigate the CP violation effects in the long-baseline neutrino 
oscillation experiments in the four-neutrino model described in Sect. II. 

We adopt the most general parametrization of the mixing matrix U for Ma- 



jorana neutrinos 1 20], which includes six mixing angles and six phases. The ex- 
pression of the matrix is too complicated to write it down here. So, we cite 
only the matrix elements which are useful for the following numerical analyses; 

Uel = C01C02C03, ?7e2 = Co2Co3S^oH ^e3 = Co3S^q2, ?7e4 = -5^03, t/^3 = —S^02-^d03S'^l3 + 
— Co3'5^i3, l/t3 — Ci3S^Q2Srf03'S^23 Co2S^;^2"^rfl3'S^23 "I" '^02Cl2C23, Slid 
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Ur4: = co3Ci3sjj23, wheie Cij = COS 9ij and Sdij = Sije^^^^ = sin%e"^'J |]2y], and 
^01, ^02, 6*03, 6112, 6*13, 6*23 are the six angles and 5oi, ^02, ^03, 5i2, ^13, ^23 are the six 
phases. As stated in Sect. II, three of the six oscillation probability differences 
are independent so that only three of the six phases are determined by the mea- 
surements of the CP violation effect in the neutrino oscillations. In this sense, our 



analyses apply both to the Dirac and Majorana neutrinos ||22||. 

On the basis of this parametrization, we obtain the constraints on the mixing 
angles and phases by using the constraints on the mixing matrix elements derived 
in the previous section. First, the constraint of Eq.(|iy) leads to 

4^02 + ^03 < 0.025. (17) 
This unequality means at least So2,Sq3 < 0.025. The next constraint of Eq.(|D of 

I - so2So3Si3e-'^'°'~'''^''''^ + co2Ci3Si2e-''''\^cl^sl^ > 0.205 (18) 

leads to 

suclssl, > 0.205 (19) 

due to the smallness of S02 and S03. The third constraint of Eq.(0) gives the 
following expression, 

IC02S02C03S12C13 + 0020035035136"^' | = 0.02 - 0.16, (20) 

where 6i = 602 — ^03 — ^12 + ^is- This constraint proves not to bring any constraint 
on the phase 61, if we use Eqs.(p!7D and (p^. The fourth constraint of Eq.(|l6[) is 
expressed as 

I CQ2C12S12C13C23 — Co2So2So3'5l2C^3S23e '('^^+'^2) _ Co2So2So3Cl2'Si3C23e'^' 

+ Ci3Si3S23(Cq3 — Cq2S^2 + '5o2'^03)^ '^^1 < 0.28, (21) 

where 62 = 5i2 — ^13 + ^23- By using Eqs. ( p!7| ) and (|1^), no constraint on ^i, 
and the fact of the large angle mixing in z/^ oscillation for the atmospheric 
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neutrino anomaly which leads to the nearly maximal mixing in the angle 623, the 
constraint of Eq. ( pi]) gives no constraint to the phase 62- 

So, in summary, we derive the two constraints of Eqs. ([T7| ) and ([191) on the 
mixing angles and no constraints on the two phases of 61 and 62- 

Using these two constraints on the mixing angles, we will calculate the differ- 
ences of the oscillation probabilities between the CP-conjugate channels for the 
long-baseline neutrino oscillations. As stated before, only three of the six proba- 
bility differences among the four neutrinos are independent so that three of the 
six phases are relevant here. However, only two phases dominantly affect the dif- 
ferences as is shown by the leading terms relevant to the long-baseline oscillation, 
which are given in the following, 

AP^e ^ 4Co3Co2So2So3-Sl2Cl3Si3 siu 5i siu 2A43, 

APer ^ 4C03C02S02S03C13S23 [-C12C23 sin((5i + 62) + S12S13S23 sin Si] sin 2A43, 

AP^^ ~ 4Co3Co2Ci3Si3S23[Co2Cl2Si2Ci3C23 siu §2 + S02S03S12S23 siu 5i 

- S02S03C12S13C23 sin(5i + ^2)] sin 2A43, (22) 
APe^ ~ 4C03C02S02S03C13C23 [C12S23 sin(5i + ^2) + S12S13C23 sin 61] sin 2A43, 

AP^s - 4Co3Co2Ci3Si3C23[co2Cl2Sl2Ci3S23 siu 62 - S02S03S12C23 siu 61 

- S02S03C12S13S23 sin(5i + ^2)] sin2A43, 

APrs ~ -4C03C02C12C13C23S23 [C02S12S13 sin 62 + S02S03C13 sin(5i + ^2)] sin 2A43, 

where 61 and 62 are the linear combinations of 6i/s as stated before. We take the 
range of phases as < 6ij < 27t and the range of mixing angles as < 6ij < tt 
so that Sij's can be taken only positive and Cij S C3ill be taken both positive and 
negative. Since Eq. (IT^) means that the angles 6q2 and 6'o3 are very small and 
Eq. (^) leads to 5^28111^2^^13 > 0.82 which means that S12 is in the range of 
0.9 < S12 < 1.0 and the angle 6^13 is around 7r/4, AP^e and APe,- are expected 
from Eq.dH) to be very small and AP^,- is to be able to take a sizable magnitude. 
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Table 1: Phase 5i-dependence of the oscillation probabilities Pa/s and their 
differences APa/3 for the long-baseline experiment. P^^ = P(t'a ~^ ^p) and 
APa^ = P{va yp) - Pi^a i'p)- ^1 is in degree. 
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p 
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0° 


0.032 


0.004 


0.148 


0.000 


-0.007 


-0.269 


45° 


0.036 


0.008 


0.148 


0.015 


0.002 


-0.277 


90° 


0.032 


0.015 


0.147 


0.022 


0.011 
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0.022 


0.021 


0.147 


0.015 


0.013 


-0.282 


180° 


0.011 


0.022 


0.148 


0.000 


0.007 


-0.275 


225° 


0.006 


0.018 


0.148 


-0.015 


-0.002 


-0.267 


270° 


0.011 


0.011 


0.148 


-0.022 


-0.011 


-0.261 


315° 


0.021 


0.005 


0.148 


-0.015 


-0.013 


-0.262 



Table 2: Phase 52-dependence of the oscillation probabilities Pap and their differ- 
ences APq/3 for the long-baseline experiment. 62 is in degree. 
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0.032 
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0.022 


0.003 


-0.008 


45° 


0.032 


0.009 


0.124 


0.022 


0.006 


-0.201 
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0.032 


0.015 


0.147 


0.022 


0.011 


-0.283 


135° 


0.032 


0.023 


0.248 


0.022 


0.016 


-0.205 


180° 


0.032 


0.028 


0.367 


0.022 


0.018 


-0.014 


225° 


0.032 


0.028 


0.433 


0.022 


0.016 


0.179 


270° 


0.032 


0.022 


0.409 


0.022 


0.011 


0.261 


315° 


0.032 


0.014 


0.309 


0.022 


0.006 


0.184 



In the following, we calculate the oscillation probabilities P(z^ci —>■ z/^) and their 
differences AP^^ by using the rigorous expressions of Eqs. (|]) and (0). 

The probabilities of P(t'|j t'e) and P(t7i — > t^e) as functions of the phase 
with (^2 = 7r/2 fixed are shown in Fig.l and those of P(z/e ^r) and P{i7e ^r) 
as functions of the phase 62 with 61 = 7r/2 fixed are shown in Fig. 2 for the 
values of the parameter set of angles and phases; S02 = sqs = 0.11(co2 = C03 = 
0.994), S12 = 0.91(ci2 = 0.415), si3 = 0.67(ci3 = 0.742), soi = S23 = 1/^2(001 = 
C23 = l/"\/2) and 5oi = ^02 = (^os = S12 = 0, which are chosen so as to give the 
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probability differences as large as possible within the parameter ranges allowed 
by the constaints of Eqs. (|T^ and (T^). The magnitude of these probabilities 



is at most 0.04 as shown in Figs.l and 2. Therefore, the probability differences 
AP(z/^ Ue) and AP(z/e — » z/^) are at most ±0.02 as shown in Fig.3 for the same 
parameter values. On the other hand, Pii'^ Vt) and P{y^ ^v) can rise to 
as large as 0.40 — 0.45 as shown in Fig.4 and AP(z/^ u^) can attain as large as 
±0.28 as shown in Fig.5 for the same parameter values. These facts agree with 
the above-mentioned expectations. The angle ^23-dependence of AP(z/^ Ur) 
and AP(z/e z/^) is shown in Fig. 6, where the phases 5i and 82 are taken as 
7r/2 and the values of the other angles and phases are the same as the above. 
We display the phase 5i-dependence in Table 1 and the phase 52-dependence in 
Table 2 of P(z/^ -> z/e),P(z/e z/^),P(z/^ z/^), AP(z/^ ^ Ue),AP{ue -> Ur) and 

Here we comment on the matter effect on the oscillation probability difference 
in the four-neutrino model. By using the Minakata-Nunokawa procedure^], the 
probability difference with the matter effect is expressed for the long-baseline 
Ua Up oscillation in the four-neutrino model with mass scheme of the two 
nearly degenerate pairs separated with the order of 1 eV as follows, 

APo,i3 ~ ilmiU^^U^sU^^Ufii) cos P34 sin i j 



+ 4Re{U*^UM3Up4) sin P34 sin 



2E 



8 E ReiUUU5V)^,,, cos' {^) sin^ (^^) , (23) 



j>i 



where 



B,, = {\Uei\' - \UejY)aL ± {\UsiY - \UsjV)a'L, (24) 
+ U:,6V^,U;jUp, + 6V:,U^,U*p^Up,. (25) 
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In Eq.(^^, the quantity a represents the matter effect for z/g and we take a = 
1.04 X 10~^^ eV for the constant matter density of 2.7 g/cm^ and a' represents 
the one for and we take a' = a/2 ||15|. In Eq . (^Sf) , 6Vai is given p| by 



= E -^U^,{U:^Ue^a + U:^Us^a!). (26) 

In Eq.(^), the first term represents the genuine CP- violation effect corrected 
by the matter effect, the second term does the CP- violation effect coming from 
the phase evolution of the neutrino wave function in the matter, and the third 
one results from the corrections to the mixing matrix V due to the existence of 
matter 0. 

We estimate these matter effects for the z/^ oscillation. The first term 

of Eq.(p3D is almost the genuine CP- violation effect, since the magnitude of the 
matter effect i?34 is at most 1 x 10^'^ for the above-mentioned parameter values 
of the mixing angles and phases. The second term is approximately 0.4 x 10^^, 
since Re([/*4t/^3t/*3t/T-4) ~ 0.1 and sin ~ 1 x 10^^. This should be compared 
with the possible maximum value of the genuine CP-violation effect displayed in 
Fig. 5 and Table 1, that is, |AP(z/^ ^ z/^)| ~ 0.3. The third term of Eq.(P5|) is 
expressed as 

/A/\/f2r\ ( /\m'^j\ 

-8 MUUU5V),r;,sm' ( -8ReiUUU5V),,,usm' (^^1 , 

(27) 



since cos (Bij/2) ~ 1.0. The coefficient of the first term of Eq . (|27|) is given by 

J2 Re iUUU6V)^r;ij 

i=l,2,j=3,4 

= 4|^Re[([/;3[/.3 + u;,Ur,){iu,3u:, + u,,u:,)iu:,u,s + u:,u,,) 

+ liu,su:, + u,,u:,)iu:,Uss + k.u,,)}], (28) 

where the relation a' = a/2 is used, and the terms with l/Am^ and l/Am^^i.^^ do 
not appear due to the symmetry of the mass scheme of the four neutrinos adopted 
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in our model. The coefficient of tlie second term of Eq.(|27|) is given by 



Re {UUUSV)^r;34 

2Ea 



+ 2^33^ + 2f/;3f/Mf/.3?7;4(|f/e3r + WeA' + l^Ws^V + 2 1^341')] 

+ '^M{u,3u;,{\Ur,\^ - |t/.4r) + Ur,uu\u,,\^ - \u,,\^)}{u:^Ue, 

+ \u:,U,,)]. (29) 

The magnitude of Eg . (pSf) is estimated to be —0.9 x 10"^, and the magnitude of 
Eq.dl) is to be 1.1 x 10"^ 0.51 x lO^^ 1.3 x lO^^ for Am^ = (1.0, 2.0, 5.0) x 10-^ 
eV^, respectively. So, the third term of Eq. (p3D , that is, Eq. (|27|) for 
oscillation is again negligibly small as compared with the possible maximum value 
of the genuine CP- violation effect. So, the matter effect can be totally neglected 
in the Vy^ —* oscillation in the four-neutrino model with mass scheme of the two 
nearly degenerate pairs separated with the order of leV, as was generally studied 
for any channels in ref . p2| . 



As can be seen in Figs. 5 and 6 and in Tables 1 and 2, CP violation could 
be observed as the probability difference between the — > Uj. and Ur 
oscillations in the four- neutrino model. So, we show in Figs. 7 and 8 the oscillation 
probabilities -P(t'^ — > Vt) and P{Vf_i —>■ Vt)-, and their difference AP(z/^ jj^) as 
functions of Am?/E [eV^/GeV], respectively, for the long-baseline experiments of 
the MINOS m and CERN-ICARUS |3|] types, where E is the neutrino energy. 



In Figs. 7 and 8, we have assumed the baseline length as L = 730 km. We can 
observe from Fig. 8 that if the beam energy is taken as 7 GeV, magnitude of the 
CP violation effect for the z/^ channel could attain as large as |AP| ~ 0.22 

in the case of Am^ ~ 3.5 x 10^'^ eV^. Incidentally, we display the probabilities 
P{^fj. ~^ ^e) and P{vy — > Ve) in Fig. 9 and AP(z/^ — > z/g) by a dashed curve in 
Fig. 8. We show in Table 3 the Am^/ii^-dependence of Piv^i —>■ z/e),P(z/e 
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Table 3: Am^/£^-dependence of the oscillation probabilities Pap and their differ- 
ences APap for the long-baseline experiment of L = 730 km. Paf3 = P{i'a ^p) 
and AP^p = P{u^ up) - P{i7^ i7p). Aw?/E is in 10-=^ eVVGeV. 



Am^/E 




p 

er 


p 




APe. 


AP^, 


0.1 


0.023 


0.007 


0.068 


0.004 


0.002 


-0.051 


0.2 


0.025 


0.008 


0.053 


0.008 


0.004 


-0.101 


0.3 


0.027 


0.009 


0.046 


0.011 


0.006 


-0.147 


0.4 


0.028 


0.010 


0.047 


0.014 


0.007 


-0.188 


0.5 


0.030 


0.011 


0.056 


0.017 


0.009 


-0.223 


0.6 


0.031 


0.012 


0.072 


0.019 


0.010 


-0.251 


0.7 


0.032 


0.013 


0.096 


0.021 


0.010 


-0.270 


0.8 


0.032 


0.014 


0.126 


0.021 


0.011 


-0.281 


0.9 


0.032 


0.015 


0.161 


0.022 


0.011 


-0.282 


1.0 


0.032 


0.016 


0.201 


0.021 


0.010 


-0.274 


1.2 


0.030 


0.016 


0.287 


0.018 


0.009 


-0.231 


1.4 


0.028 


0.016 


0.374 


0.012 


0.006 


-0.158 


1.6 


0.024 


0.014 


0.450 


0.005 


0.002 


-0.064 


1.8 


0.020 


0.012 


0.504 


-0.003 


-0.001 


-0.038 


2.0 


0.016 


0.010 


0.530 


-0.010 


-0.005 


0.135 


2.5 


0.011 


0.005 


0.460 


-0.021 


-0.011 


0.279 



Vr), P{i^n i^t), AP(i/^ Ue), AP{ue ^ z^r) and AP(i/^ Ur) for L = 730 km. 

V Conclusion 

We have derived the constraints on the neutrino mixing matrix by using the 
data from the solar neutrino deficit, atmospheric neutrino anomaly, LSND oscilla- 
tion experiments, Bugey experiment and the CHORUS and NOMAD experiments 
along with the other accelerator and reactor experiments in the four-neutrino 
model with mass scheme of the two nearly degenerate pairs separated with the 
order of 1 eV. We have used the most general parametrization of the mixing ma- 
trix with six mixng angles and six phases applicable to both Majorana and Dirac 
neutrinos and have obtained the two serious constraints about the four of the six 
mixing angles and no constraints on the phases. 
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By using these constraints, we have calculated the oscillation probabihties of 
P{^n ^e), P{^e J^t) and P(t'^ — > z^r) and have investigated CP violation in the 
long-baseline neutrino oscillations of AP(i/^ i/g), AP(i/e —>■ 1/^) and AP(i/^ — > 
Vr). The quantity AP(i/^ — > Vt) is found to be able to attain a value as large as 
±0.28 due to the large mixing between and v^. and the mass scheme of the four 
neutrinos and, therefore, it could be observed in the long-baseline experiments. 
We have shown that the contribution to AP(i/^ — > i/^) from the matter effect is 
at most 0.01, 0.005 and 0.001 in magnitude for Am^ = 1.0 x 10-^ 2.0 x 10"^ and 
5.0 X 10~^eV^, respectively. So, we can conclude that the matter effect is neghgibly 
small in comparison with the possible maximum value of the genuine CP- violation 
effect of |AP(z/^ ^^t)\ = 0.28. 
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Figure captions 



Fig.l. The oscillation probability of Vfj, — > i^e (solid curve) and — > z7e (dashed 
curve) with respect to the phase 5i of the mixing matrix for the long-baseline 

experiment. The other angles and phases arc fixed as So2 = Soa = 0.11(co2 = C03 = 
0.994), S12 = 0.91(ci2 = 0.415), si3 = 0.67(ci3 = 0.742), sqi = S23 = l/V2(coi = 
C23 = 1/v^), 5oi = S02 = 5o3 = 5i2 = and 82 = 7r/2. 

Fig. 2. The oscillation probability of — > i^r (solid curve) and z/g — > (dashed 
curve) with respect to the phase 82 of the mixing matrix for the long-baseline 
experiment. The other angles and phases are the same as in Fig.l except for 
5i = 7r/2 fixed. 

Fig.3. The probabihty difference AP(i/^ — > Ue) (sohd curve) and AP{ue 
i^r) (dashed curve) with respect to the phase Si for the long-baseline experiment. 
The other angles and phases are the same as in Fig.l. 

Fig. 4. The oscillation probability of Uj. (solid curve) and 1/^ z7^ (dashed 

curve) with respect to the phase S2 for the long-baseline experiment. The other 
angles and phases are the same as in Fig. 2. 

Fig.5. The probability difference AP(i/^ — > u-,-) (solid curve) and AP(i/e — > Ur) 
(dashed curve) with respect to the phase S2 for the long-basehne experiment. The 
other angles and phases are the same as in Fig. 2. 

Fig. 6. The probability difference AP(z/^ v^) (solid curve) and AP(z/e — t-'t) 
(dashed curve) with respect to the angle ^23 for the long-baseline experiment. The 
other angles and phases are the same as in Fig.l except for 5i = 7r/2 fixed. 
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Fig. 7. The oscillation probability of u^j, — > Ur (solid curve) and 17^ — > (dashed 
curve) with respect to Arri^/E for the long-baseline experiment with the distance 
of L = 730 km. The angles and phases are the same as in Fig.l except for 5i = 7r/2 
fixed. 

Fig.8. The probability difference AP(i/^ — > i/^) (solid curve) and AP(i/^ — > 
fg) (dashed curve) with respect to Arrr^/E for the long-baseline experiment with 
the distance of L = 730 km. The angles and phases are the same as in Fig. 7. 

Fig.9. The oscillation probability of z/^ z/g (solid curve) and z/^ (dashed 

curve) with respect to Arv? j E for the long-baseline experiment with the distance 
of L = 730 km. The angles and phases are the same as in Fig. 7. 
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